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a 

The free photon dispersion relation is a reference quantity for high precision tests of Lorentz Invariance. 
We first outline theoretical approaches to a conceivable Lorentz Invariance Violation (LIV). Next we ad- 
dress phenomenological tests based on the propagation of cosmic rays, in particular in Gamma Ray Bursts 
(GRBs). As a specific concept, which could imply LIV, we then focus on field theory in a non-commutative 
(NC) space, and we present non-perturbative results for the dispersion relation of the NC photon. 
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1 Lorentz Invariance Violation (LIV) 



C*~) ■ Lorentz Invariance (LI) plays a central role in relativity: it holds as a global symmetry in Special Relativity 

Theory, and as a local symmetry in General Relativity Theory. We refer to the former in order to be able to 
describe particle physics in terms of quantum field theory, where a field $ (scalar, vector, tensor or spinor) 
is transformed globally in some representation D of the Lorentz group SO(l, 3), 



$0) -> U(A)^{x)U(A) = D(A)$(A _1 x) , U unitary , A G 50(1,3) . (1) 



If a local field theory obeys LI, then it is also CPT invariant JT]. If LI is broken, CPT may be intact or not. If, 
however, CPT is broken, then a LI Violation (LIV) is inevitable |2j. Therefore LI can be probed indirectly 
by testing CPT predictions, such as the equivalence of the mass and magnetic moment of particles and 
their antip articles. In particular the masses of K° and K° coincide to a relative precision < 8 • 10~ 19 |3). 

Specific direct LI tests are even more stringent. In the framework of the Standard Model Extension 
(SME) 21, local terms are added to the Standard Model Lagrangian, which can cause a spontaneous LIV, 
e.g. 

C = i^d^ - g^cfnp - ^(g'G^ + g"H^ 5 )^d v ^ + ... (2) 

Spontaneous Symmetry Breaking (SSB) of the Higgs field <fi G C 2 leads to the fermion mass. In analogy 
the new tensor fields G and H may also undergo SSB, which affects LI (unlike the SSB of the Higgs field). 
For a simple change of the observers frame, G and H are transformed as well, so that C remains invariant; 
but the fermion -0 perceives them as an invariant background. LIV is then manifest for instance through 
a modified fermion dispersion relation, depending on the Yukawa-type couplings g' , g". Kostelecky et 
al. identified more than 100 parameters of this kind (for a review, see Ref. J3]), which are claimed to 
be compatible with the usual properties of the Standard Model except for LI (and in part CPT) [6]. This 
corresponds to a low energy expansion, after lifting the usual LI condition. This collaboration interprets 
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the photon as a Nambu-Goldstone boson Q. All LIV parameters emerging in this framework are experi- 
mentally bounded by O(10 -27 ) (in a dimensionless form, which factors out the Planck scale) J8). If some 
of these parameters are non-zero, they still allow for sizable LIV effects at tremendous energies, though 
this scenario poses a severe hierarchy problem J5J- 

In an effective approach, Coleman and Glashow considered Maximal Attainable Velocities (MAVs, 
attained in the high energy limit), which could slightly differ from the speed of light, depending on the 
particle type ifTOl . They were most interested in the impact on cosmic rays (which attain the highest 
energies in the Universe). In particular the prediction that proton rays cannot exceed the Greisen-Zatsepin- 
Kuz'min (GZK) energy of Eqzk ~ 6 • 10 19 eV over a distance of more than ~ 100 Mpc in the Cosmic 
Microwave Background (CMB) (mainly due to photopion production via a A-resonance) could be evaded 
by tiny MAV differences. However, recent observations favour the validity of the GZK energy cutoff IfTTI : 
up to Eqzk the cosmic ray flux follows closely a behaviour oc E~ 2J , but beyond Eqzk there is an extra 
suppression compared to this power law. This suggests that LI holds even under boosts of Lorentz factors 
like 7 = Eqzk /mpioton ~ 10 11 , i.e. 6 orders of magnitude beyond the 7-factors accessible in laboratories, 
which further constrains the MAV differences lfT2l . 

This issue is reviewed in Ref. (T3), covering theoretical and phenomenological aspects and the current 
status, along with MAV effects in cosmic rays with respect to decays, Cerenkov radiation and neutrino 
oscillation [;14j (these effects were addressed already in Ref. liTOI . and reviewed earlier in Ref. 03)). 



2 Cosmic 7-rays 

The highest energies in cosmic 7-rays have been observed around E 1 ~ 50 TeV. They were emitted from 
the Crab Nebula lfl6l (a remnant of a supernova of the year 1054, at a distance of 2 kpc). With certain 
assumptions about the origin of these rays (synchrotron radiation of high energy electrons and positrons, 
followed by powerful inverse Compton scattering) one obtains stringent LIV constraints lfT31 . 

The most energetic 7-rays reaching us from outside our galaxy originate from blazars, a subset of 
the Active Galactic Nuclei. Two prominent blazars are denoted as Markarian 501 and 421 at distances 
of about 157 Mpc and 110 Mpc, which had emitted photons that arrived at Earth with E 1 ~ 20 TcV 
IfTTI resp. E 1 ~ 10 TeV IfTSl . These multi-TeV 7-rays have been considered a puzzle similar to the 
GZK cutoff (which seemed to disagree with the data in the 20 th century): they give rise to the creation 
of electron/positron pairs in the CMB , which can perform inverse Compton scattering on further CMB 
photons. This leads to a cascade of high energy photons, which are attenuated in this way. The penetration 
depth of these cascades in the CMB, above certain energy thresholds, is reviewed for instance in Ref. ||T9l : 
it is sensitive to details of the extragalactic magnetic fields and the radio background, which are not known 
precisely. In analogy to the GZK energy cutoff, also here LIV effects — and in particular MAV differences 
— have been proposed as an explanation for the possibly puzzling 7-radiation from blazars l20l . However, 
also in this case it is not clear if they really represent a puzzle — e.g. the analysis in Refs. ||2TI suggest that 
their understanding does not require new physics. 

Since 1967 we also know about the existence of Gamma Ray Bursts (GRBs), which are emitted in 
powerful eruptions. Temporarily they yield the brightest spots in the sky. Typically they last for a few 
seconds or minutes, and they involve photons in the energy range E 1 ps (10 4 . . . 10 8 ) eV. About their 
origin we only have speculations, such as the merger of black holes or neutron stars, along with fireball 
shock models for the GRB generation. In any case their existence is suggestive for a high precision test 
of the photon dispersion relation, which represents a purely kinematic facet of LI (22 1- A deviation from 
the linear dispersion relation, i.e. a photon group velocity v 1 (E^) ^ const. , can be obtained most simply 
by introducing a photon mass m 7 > (through a modified Higgs mechanism, keeping the photon MAV 
unaltered). However, the bound based on laboratory experiments, m 7 < 10~ 18 eV Q, is so tiny that its 
impact would still not be detectable in GRBs 11221 . Hence this approach is not attractive in this context and 
we return to m~ = 0. 
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Instead we switch to the simple effective ansatz (with the photon speed c = 1 attained at low E y ) 



I - El 

M ' 



(3) 



where M is a very heavy mass parameter (which could emerge from some unspecified "quantum gravity 
foam"), and we consider only 0{E 1 /M). Ref. Il23l confronted this ansatz with the data of 35 GRBs 
measured from three satellites, and observed that photons with higher energy tend to arrive earlier. They 
parametrised the observed time delay At between energy channels, which differ by AE, as 



At 

T+z 



z+a-LW K(z) , a L iv 



AE 
H M 



Hq : Hubble parameter, z : redshift . (4) 



The parameter d SOU rcc captures a possible relative delay already at the source, and K(z) is a function of the 
redshift, which is derived in Ref. 12311 . The question is now whether or not the data are compatible with a 
vanishing LIV parameter. Fig.[T](on the left) shows that they do certainly not rule out oliv = 0, although 
the best linear fit has a slight slope. The plot on the right specifies the windows for the parameters G? SOU rcc 
and aLiv with 68 % resp. 95 % confidence level (C.L.). The observation that a^iv cannot deviate much 
from zero implies 



\M\ > 1.4 



10 25 eV 



0.001 M Planck (with 95 % C.L.) . (5) 
This result appears robust due to the large data set considered. Similar studies based on single GRBs or 
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Fig. 1 Evaluation of the parameters in eq. {4]l from GRB data of three satellites. The plot on the left shows the left- 
hand-side of eq. 10 against the function K(z). The plot on the right illustrates the resulting bounds for the parameters 
dsource (relative delay at the source) and <xliv (LIV effects) with 68 % resp. 95 % C.L. (plots adapted from Ref. [23 1). 



blazar flares (again involving Markarian 501) did not report evidence for LIV either, and derived bounds 
for \M\ even up to 0(0.01) • M P i anck |24| . 

As another purely kinematic effect — which does not require any assumptions about the interactions in 
a LIV modified form of QED — there could be vacuum birefringence, if w 7 is slightly helicity dependent, 

k 

v^(k) = 1 ± (k = \k\ : momentum, Afh : heavy mass, ± : hclicitics) . (6) 

In this case linearly polarised light should be depolarised completely after a long enough time of flight. 
However, cosmic 7-rays from quasars at a distance of 300 Mpc, with wave lengths of (400 . . . 800) nm 
were observed, essentially with a linear polarisation (less than 10° deviation from the polarisation plane) 
1251 . This implies for that LIV scenario even a strongly trans-Planckian bound: \M^\ > 10 4 Mpi anc k l26l . 

Still better precision is expected soon from GRB and blazar flare photons with E 1 sa 8 keV . . . 300 GeV, 
to be detected by the Fermi Gamma-ray Space Telescope 1271 . which arrived in orbit in June 2008. 
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In addition the Chinese-French Space-based multi-band astronomical Variable Objects Monitor mission 
(SVOM) (28) plans to monitor about 80 GRBs a year, starting in 2013. 

Also the subjects of this section are reviewed more extensively in Refs. II131I15I . 

3 £7(1) gauge theory in a non-commutative (NC) plane 

We now address another theoretical concept for LIV, which is fundamentally different from the low energy 
effective approaches (4][10) that we referred to in Section 1 . We first consider a Euclidean plane given by 
Hermitian coordinate operators xi, £2, which do not commute. This yields a purely spatial uncertainty 
relation, 



ip,x„] = iOe^u , Axi • Ax 2 > -\9\ . (7) 



£1 

Here we assume 9 = const, (which can be generalised, of course). \Z\6\ now represents a kind of "minimal 
length", as a new constant of Nature, like the parameters 1/M or in Section 2. It might have an 

interpretation as the event horizon of a mini black hole in a (hypothetical) attempt to measure distances 
of 0(Lpi ari ck) (which would require a tremendous energy density) l29l . Field theory on such a space is 
necessarily non-local over this range. This gives rise to the mixing between UV and IR divergences, which 
makes perturbation theory mysterious (for a review, see e.g. Ref. BOl and references therein). 

Here we refer to a fully non-perturbative approach, which starts by introducing a lattice structure. Let 
the spectra of x\ and £2 be discrete multiples of some "lattice spacing" a (although there are no sharp 
lattice sites, since x\, £2 cannot be diagonalised simultaneously). On the other hand, the momentum 
components k\ , k 2 commute, and they obey the usual periodicity over the Brillouin zone, 

e ifc M i M _ e i(fc M +27r/a)£ M . ^ _ g i(fc t , +2ir/a)£ M -ifc„i„ _ j^ i7T0(fc 2 -fci )/a u 

2a M 

For fixed parameters 9 and a the momenta turn out to be discrete, i.e. the lattice must be periodic. Let us 
assume periodicity over a N x N lattice, hence the fc M are integer multiples of 2n/(aN), and from the 
computation above we infer 

9 = -Na 2 . (8) 

In our studies we were therefore interested in the Double Scaling Limit (DSL), which takes a — > and 
N — > 00, while keeping Na 2 = const. This simultaneous UV and IR limit leads to a continuous NC plane 
of infinite extent. The necessity to link these limits is related to the UV/IR mixing mentioned above. In 
other regularisation schemes it is less obvious how to achieve this delicate balance. 

It is equivalent to return to coordinates x^, but perform all field multiplications by a star product, e.g. 

4>{x) * tp{x) := 4>{x) exp Q \ 9 e^Sj ip(x) . (9) 
Thus the NC U(l) pure gauge action takes the star-gauge invariant form 

S[A] = \j d2x F P-y * F nv 1 F ^ = d P- A » ^ d » A n + [ 9( A ^ *A V -A V * A„) , (10) 

which involves a Yang-Mills term. Its lattice discretisation does not appear simulation-friendly; note that 
compact link variables would have to be star-unitary. However, there is an exact map onto a matrix model 
in one point — the twisted Eguchi-Kawai model PTIl — which is Morita equivalent [32]]. This renders 
simulations feasible, and it also provides a sound definition of Wilson loops. These Wilson loops are 
complex in general, but the lattice action is real positive, since it sums over both plaquette orientations. 
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Our Monte Carlo simulations — using a heat-bath algorithm — revealed the following features 031 : at 
small area the Wilson loops (W) of all shapes are real and they follow an area law behaviour. At large area, 
however, a complex phase sets in, which obeys for many shapes the simple relation arg(W) = arca/6>. If 
we formally identify 1/9 with a magnetic field across the plane, then this is the behaviour of the Aharonov- 
Bohm effect. This identification can be conjectured on tree level, and it has been used extensively in the 
literature; for instance, it is inherent in the Seiberg-Witten map [34|. In the case discussed above, however, 
this relation only shows up as a non-perturbative effect. 

Moreover the Wilson loops of large area were observed to be shape dependent, so the symmetry under 
Area Preserving Diffeomorphisms is broken, including the subgroup of linear unimodular transformations, 
SL(2, R) 051 . That symmetry was a corner stone in the solution of Yang-Mills theories in a commutative 
plane. We conclude that the same technique cannot be applied in the NC plane, so an analytical solution 
would be difficult, but the system may have a rich structure. 

4 The fate of the non-commutative photon 

It is straightforward to extend action ( TTOb to higher dimensions and any U (n) gauge group in the funda- 
mental representation, but the formulation of SU(N) gauge theories in NC spaces fails. Therefore 4d NC 
U(l) gauge theory is appropriate to investigate a link to particle phenomenology. 

On the perturbative side, the 1-loop calculation of the effective potential suggests a photon dispersion 
relation of the form 0"6l 

E 2 =p 2 +C-^— , pop := -9""?,,^% > , (C^const., = -i[x^,x v ]), (11) 

pop F 

which contains an IR divergence at finite 0. Moreover we see that the transition to a vanishing non- 
commutativity tensor, <d — >• 0, is not smooth (which invalidates an expansion in small ||0||). Therefore 
one could even feel tempted to claim a lower bound for |jO|| — '/it is non-zero — from the fact that no 
deviation from the linear dispersion has been observed down to small photon momenta 071 . 

However, the situation becomes even more dramatic due to the observation that the above 1-loop cor- 
rection is actually negative, C = —2/ir 2 [38]Q This looks worrisome indeed: the theory seems to be IR 
unstable. Without a ground state the system does hardly have a physical interpretation. The negative IR 
divergence is avoided in a supersymmetric version, but here we are interested in the NC photon without a 
SUSY partner. We addressed the issue of its stability in a fully non-perturbative study [40]. 

To this end we considered a Euclidean space consisting of a commutative and a NC plane. The for- 
mer contains the Euclidean time. This is essential in order to preserve reflection positivity as one of the 
Osterwalder-Schrader axioms, which are required for the transition to Minkowski signature. Once time is 
commutative, non-commutativity can be restricted to one plane by means of rotations. This setting is the 
safest NC formulation beyond perturbation theory, and it can be viewed as minimal non-commutativity; 
there is again just one non-commutativity parameter 9, as in eq. 10. 

We regularise the commutative plane by a standard LxL lattice, and the NC plane with an A x A lattice 
(L « A), which is again mapped onto A x A matrices in the twisted Eguchi-Kawai model. The first goal in 
our numerical study [40] was the identification of a DSL. The results for a set of observables (Wilson loops 
and lines and their 2-point functions, matrix spectra as coordinates of a dynamically generated geometry) 
stabilised to a good approximation if different A and (3 oc l/g 2 were tuned as 

N/f3 2 = const. aocl//3 and 9 oc N/0 2 . (12) 

As an example we show in Fig.|2]the real part of square shaped Wilson loops in different planes, at a fixed 
ratio N/ (3 2 = 20. In the commutative plane and in the mixed planes (W) is real due to reflection symmetry 
on the commutative axes. In the NC plane we see an oscillating behaviour of the real part at large area, 
which is due to a rotating complex phase, similar to the 2d case f33l that we described in Section 3. 

' On the other hand, a positive IR divergence has been observed numerically in the NC A0 4 model 1391 . 
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Fig. 2 Wilson loops measured in NC QED4 at N / /3 2 = 20 in the commutative plane, a mixed plane and the NC plane 
(from left to right). In the first two planes (W) is real, whereas the oscillation of Re (W) in the NC plane is due to a 
rotating complex phase (similar to NC QED2 133V ). 




P P 

Fig. 3 The expectation values of open Wilson lines of length 2 and 4 lattice spacings ( on the left/right) for TV = 15, 25 
and 35 (curves connecting bullets, triangles and squares, with hysteresis at the transition to the weak coupling phase). 
At intermediate gauge coupling we discovered a phase of spontaneously broken translation symmetry. 



The next goal was the exploration of the phase diagram. For this purpose we measured expectation 
values of open Wilson lines P^ (n) in the NC plane, with a length Q^Pv = nap, of a few lattice spacings 
(\fi\ = 1, n G IM). The P^n) are gauge invariant in NC space; note that star-gauge transformations 
are non-local themselves. Since these open Wilson lines carry finite momentum, they are suitable order 
parameters to detect the spontaneous breaking of translation symmetry. Results for N = 15, 25 and 35 are 
shown in Fig. [3] Translation symmetry is intact at strong and at weak gauge coupling, in agreement with the 
corresponding expansions, but it is broken at intermediate coupling. This effect is purely non-perturbative; 
it was not visible in analytical studies. The transition strong/moderate coupling occurs at f3 pa 0.35 (for 
any TV), whereas the transition moderate/weak coupling roughly follows a behaviour j3 oc TV 2 (that phase 
transition is of first order, as we see in Fig. [3] from the hysteresis). The existence of a broken symmetry 
phase at moderate coupling was also confirmed in simulations of the closely related 4d twisted Eguchi- 
Kawai model 1:411. In our case we obtain the phase diagram in Fig. [4] The DSL is attained by following 
lines of constant 9. Relation ( fl~2l implies that they correspond to (3 oc \/~N, hence the DSL always leads to 
the broken phase. Since we found stability for a number of observables as we approach this limit, it seems 
to provide IR stability for the NC photon. 

At last we consider non-perturbative results for the NC photon dispersion relation. The plot in Fig. 
|5] on the left shows the behaviour in the symmetric phase at weak coupling: it is fully consistent with 
the perturbative result of a negative IR divergence. However, we saw that the DSL leads to the broken 
phase, which is therefore physically relevant. For that phase the plot on the right of Fig. [5] reveals an IR 
stable behaviour, at least for the energy depending on the momentum component in the commutative space 
direction only, E(p = P3) Pi=P2= q. Therefore the photon could survive the transition to an NC world 
after all. Then one could again try to interpret it as the Nambu-Goldstone boson of the SSB of Poincare 
symmetry, similar to the interpretation in the context of the SME Q (cf. Section 1). 
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Fig. 4 The phase diagram for QED4 with a NC plane (in lattice units): between the strong coupling phase (f3 < 0.35) 
and the weak coupling phase we recognise a phase of broken translation invariance. The weak/broken transition rises 



like j3 oc N , whereas the DSL for a fixed NC parameter 6 follows /3 oc V iV, so any DSL ends up in the broken phase. 
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Fig. 5 Dispersion relations E{p) for the photon in a NC space, in the weak coupling phase (on the left), and in the 
broken phase (on the right). The former result is consistent with the perturbatively predicted IR instability ( "tachyonic" 
behaviour). In the physical phase, however, the photon is massless again, corresponding to a Nambu-Goldstone boson 
of the broken translation symmetry. 

5 Conclusions 



Cosmic 7-rays enable LI tests to extreme precision. So far LIV has not been detected anywhere — we 
reviewed failed attempts based on GRB data. 

As theoretical frameworks leading to LIV we discussed on one side low energy effective theories ||4}j7| 
and an effective ansatz which introduces particle specific Maximal Attainable Velocities iflOl . Another 
approach is based on NC geometry — the presence of IR divergences in quantum field theory on such 
spaces clarifies its distinction from low energy effective theories. 

In Section 4 we were concerned with the kinematics of a NC photon. Perturbation theory to one loop 
suggests its IR instability J38), which looks like a disaster for NC field theory. The picture changes, 
however, in light of a non-perturbative consideration. 

Our numerical study revealed a phase of spontaneously broken Poincare symmetry, which is not vis- 
ible to perturbative calculations. In fact the limit to the continuum and to infinite volume at fixed non- 
commutativity (the Double Scaling Limit) leads to this phase of broken symmetry,and not to the sym- 
metric weak coupling phase. In the broken symmetry phase we observed stability for a number of observ- 
ables BUI . In particular the photon dispersion relation in a commutative sub-space appears linear. However, 
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its behaviour at non-zero momentum components in the NC directions is not explored yet. In that respect, 
non-perturbative results could be confronted with phenomenological data — in particular from GRBs — 
to establish a robust bound on ||0|| in Nature. On the phenomenological side, more stringent data are 
expected from new large-scale GRB observations, for instance in the projects discussed in Refs. 12711281 . 

Our non-perturbative results so far imply that the photon may survive in a NC world — despite the IR 
disaster on the 1-loop level — without the requirement to proceed to SUSY. 
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